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TECHNICAL NOTE 3886 

AVERAGE PROPERTIES OF COMPRESSIBLE LAMINAR BOUNDARY LAYER 
ON FLAT PLATE WITH UNSTEADY FLIGHT VELOCITY 
By Franklin K. Moore and Simon Ostrach 


SUMMARY 

The time-average characteristics of Boundary layers over a flat 
plate in nearly quasi-steady flow are determined. The plate may he ei- 
ther insulated or isothermal. The time averages were found without spec- 
ifying the plate velocity explicitly except that it is positive and has 
an average value. 

Each time average Involves two groups of terms to the order consid- 
ered in the report , a time average of quasi-steady terms, and terms re- 
lated to the reduced frequency of the plate velocity fluctuations . 

The quasi-steady terms differ from the values for steady flow at the 
corresponding average velocity. These differences are reinforced By the 
frequency dependent averages for adiaBstic wall- temperature and heat- 
transfer rate. The effects oppose one another In the case of skin 
friction. 

The special case of harmonic velocity variation Is considered, and 
it is found that large amplitudes accentuate the importance of the 
frequency-dependent terms . 

Oscillating the wall to increase the heat-transfer rate is not ad- 
vantageous if the power to oscillate the plate is accounted for, unless 
the Prandtl number Is in a certain range. 


INTRODUCTION 

In many current problems of aerodynamics, unsteady motions of a sur- 
face are important . The accelerating and decelerating phases of missile 
flight and the intermittent flow in an engine during unstable combustion 
are examples. The nature of the boundary layer of these unsteady flows 
may be studied with a view to determining friction drag, surface tempera- 
ture, and rate of heat transfer through the surface. Usually, the bound- 
ary layer is so thin that It responds almost instantly to temporal changes 
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in flow conditions. Thus, the time history of the "boundary layer is a 
succession of steady states, and such a boundary layer is called quasi- 
steady. If the motion involves accelerations which are particularly rap- 
id, the quasi-steady description may require correction. 

In any case, it is often desired to assess the average effect of 
fluctuations in flow conditions. For example, if an insulated body is in 
motion with a speed which varies in time about some average value, the av- 
erage friction drag and average surface temperature may differ from values 
appropriate to steady motion at the average speed. Another example is fur- 
nished by the speculation that the net performance of a heat exchanger 
could be altered by imparting an unsteady motion to the wall. 

In the present study, as an idealized special case of the foregoing 
type of problem, a semi -infinite flat plate is assumed to be in motion 
parallel to its surface and normal to its leading edge with a flight veloc- 
ity u(t) that is always in the same direction but has a magnitude that 
fluctuates with time. The resulting boundary layer is assumed laminar and 
compressible, and the surface is either insulated or at constant tempera- 
ture. The assumed Prandtl number is 0.72, which is appropriate for air 
under normal conditions . 


The boundary layer is assumed to be nearly quasi-steady. The basic 
boundary- layer analysis is already available ’for thiB problem in references 
1 and 2, which treat the insulated plate atid constant plate temperature 
cases, respectively. The velocity profile in the boundary layer is found 
in the form 


§ - \ (a) + Z 0 f' 0 (o) + ^fi(a) + . . . + + • • D (1) 

where a. is the usual Blasius variable signifying parabolic similarity 
in the boundary layer (a complete definition will appear in a subsequent 
section) . The function F (a) is the Blasius function for steady motion 
of a flat plate. The parameters £ n govern deviations from quasi- 

steadiness . 


^n-1 


X n d°U(T) 

un+1 arjin 


. n “ > • * 


( 2 ) 


In cases of flight at substantial speed, the factor usually ensures 

that £n-l a rather small, quantity. The functions 'F(cr), fQ(cr), and 

f-^(c) are available in reference 1. A full list of symbol notation is pro- 
vided in appendix; A. 
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If the plate surface is insulated, the temperature profile may he 
written as 


^l + ^M^T)^, £ 0 ,-£i> • • •> (3a) 

where the "recovery factor" 

r * R(ff) + i Q r Q {a) + ^(cr) + . . . + ^ 0Q (a) + . . . (3b) 


and the functions R, Tq, and r^_ are tabulated in reference 1. If the 

plate is at a constant temperature at which heat transfer takes place, 
reference 2 gives 


0 - 0 . 


v 


q 2 ^ " 

riz - + 2e ~ (a^ - - ~ e ~ ) s < ff> + SopoM + S ^(e" - ej e 0 M 


. h i 


+ k (cr) + 


U 


2c P ^w “ 


3i(ff)j + . - . 


+ £oo(‘> + 2 ^Xe~Tery s oo^ ff j] + ■ ■ • 

and provides the functions H, S, h^, Sq, h-^, and s-^. 


(4) 


A study of the time-average properties of the flow represented by 
equations (l) to (4), when the speed of the plate- 1 - U fluctuates with 
time is described herein. The functions f^, r oo J ^00' an ^ s oo are 

required for this purpose, and, not being available in references 1 and 
2, are determined herein. 


A previous report (ref. 3) studies the average rate of heat transfer 
from an oscillating flat plate, but differs from the present study in that 
the plate is doubly infinite and there is no net motion through the sur- 
rounding fluid; on the average, the plate is at rest. Appropriate com- 
parisons are made herein between the results of the present study and those 
of reference 3. 


uh. the present problem, if compressibility is important, unsteadi- 
ness is restricted to enter only through motion of the plate. Thus, in a 
corresponding wind-tunnel test the model position would be varied mechan- 
ically and the tunnel flow would be held constant. The analyses of ref- 
erences 1 and 2 for a compressible fluid would not apply for a fixed model 
in a pulsing flow. 
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DERIVATION OF FORMULAS FOR AVERAGE PROPERTIES OF BOUNDARY LAYERS 

The stream function ijf for nearly quasi-steady flows is defined in 
references 1 and 2 as 

t = f(cr, t Q , £ 2 , • • .) (5) 

where 

a “ 2 ^ ^ 

The coordinates X and Y measured in a system with its origin fixed at 
the leading edge of the plate are related to those (x,y) in a coordinate 
system which is stationary in the fluid hy (see sketches (a) and (b)) 

(V) 
( 8 ) 

Equation (8) is employed to make the momentum equation Independent of the 
energy equation. The velocity in the X,Y-system is related to that In 
the x,y- system by 

u(X,Y,T) = u & (x,y,t) + U(t) (9) 





X 


Sketch (b) Coordinates fixed, in plate. 


! oyay 
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The constant C is the proportionality factor in the assumed 
viscosity-temperature variation 

4 


u =* JL = PoftVo p (10) 

y oo ^oo 

This constant may be evaluated by matching equation (10) with the Suther- 
land formula at some appropriate point, for example, at the wall. Thus 



The unsteady boundary- layer characteristics can be determined from 
equations (l) to (6) . 


Form of Unsteady Boundary-Layer Characteristics 


Insulated plate . - The unsteady boundary- layer characteristics from 
■which the average properties will be determined are presented in this and 
the subsequent section in a more extended form than in references 1 and 

2. The wall shear stress T w = {[i may be obtained in the following 

dimensionless form from equations (l), (8) , and (10), and from the state 
equation for constant pressure (p0 = Constant): 


Cf = i 


[ f " ,(0) + ? o f o(°) + V£<°> + c§ r 3o«>) 


The displacement thickness 8* may be defined^ as 


■J 

( 12 ) 

(13) 


Upon application of equation (l), the state equation (p0 =» Constant), and 
equations (3) into (13) there results 


O 

The steady-flow definition of displacement thickness is adopted here 
in. Actually, the displacement effect of the boundary layer is not proper 
ly represented by this definition if the flow is unsteady. However, an in 
vestigation to be reported elsewhere shows that the steady expression is 
part of the correct one, and since this expression yields the correct 
quasi-steady result its calculation is therefore warranted. 
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i‘ = ^|u. (20 - F) - [c 0 f 0 (“) + qfjM * 5o f OO<") ■*■■■'] 

+ m [r s ai + e °Jo* r ° d ° + c i x 40 + s °jo‘ r °° d ° + 


where 


m = (r - 1) M“ 


I 

( 14 ) 

(15) 


Isothermal plate . - The dimensionless wall-shear-stress or local- 
skin-friction coefficients for the case of heat transfer at the plate is 
identical in form with its counterpart for an insulated plate (eq. (12)). 
The constant C is evaluated "by assuming that 0 W in equation (ll) is 

the adiabatic wall temperature for the insulated plate and the maintained 
isothermal temperature in the heat-transfer case. 


The displacement thickness for this case is obtained from equation 
(l ) y the state equation, and equations (4) and (13) and is given by 


= J^/lim (2a - F) - Lf 0 (~) + S^iC-) + ^o f OO^ ) + * * 3 
l C-*» 

+ m^J'°°S da + £ 0 b q da + S-^jT s L da + ^O^J^° S 00 dcT + . • J 
+ da + C 0 j[“ h o dcJ + h l da + h °° d ° + ' * ]} 


where 


§■ 


6,, " 0 

W oo 


(IV) 


The local rate of heat transfer is given by 



( 18 ) 
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Using the definition of the Prandtl number, the state equation, and equa- 
tions (4), (8), and (10) in equation (18) yields 


<1 


- i (V S'(0) +50^(0) +§ B 6 (0) + 

+ S {C°) + • • *J + + s 6o(°) + • • )) 


(19) 


Time Averages for Arbitrary Velocity Fluctuations 


It is now assumed that the velocity fluctuates in a periodic but 
otherwise arbitrary manner, that is, 

U(t) = U m g(ost) (20) 


where <n Is the frequency of the fluctuations and g Is an arbitrary 
positive function so that 


g 



g(x) 


dT 


1 


( 21 ) 


where 

x = cut (22) 

Substituting equations (20) and (22) Into equation (2) yields 




£l 


o 2 g-.fe> 

* 3 

g 


with the frequency parameter 


given by 


S - u 
m 


(23) 


(24) 


Therefore, for arbitrary fluctuations of the velocity, the average local 

skin friction coefficient to order 2 2 is obtained by applying equations 
(20) and (22) to (24) to equation (12) and Integrating as in equation (2l). 
Since for a periodic function, 


gV = 0 


(25) 
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and also 


S 


-n g .. = ng-U+ 1 )(gi)_2_ 


(26) 


the following equation results: 

> 

2k 


^ +ffi2 [ f °° (0) + i f i <0) ] g " S/2(g,)2 } 

m (27) 

According to the remarks in the previous section, equation (27) holds for 
hoth the insulated and isothermal plates where only the value of C 
differs . 

In a similar manner the time averages of equations (14), (16), and 
(19) are 

lim (2a - F) + # W* J^R da 

{[! aa+J^T m da) g-6/2(g') 2 | 


+ S2 


for the displacement thickness of the insulated plate; 

&* = fg" 1 ^ 2 lim (2a - F) + m*g 3 / 2 S da + 2$g~ 1 / 2 H da 

t s^-[| f 1 (°°) + f 0 o (oD |] g" 9 ^ 2 (g') 2 

+ s ! da + ^ V " s 00 da ) S _5 / 2 (g') Z 

+ »4JT- h x da + J“ h ro to) g- 9 / 2 (g') 2 | U 

for the displacement thickness of the isothermal plate; and 
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for the average local heat-transfer rate of an Isothermal plate where 


and 


m = m*g 2 


(u y 

m* = m 


CpSc 


(31) 


The average value of the adiabatic wall temperature (i.e., for the 
insulated plate) is also obtained in a similar manner from equation (3) 
and is given by 


fly 

0 ® 


= 1 + 


m? 




R(0) g 2 + , 2 [r l( 0) + r 00 (0)] g^g') 2 ! (32) 


SECOND-ORDER SOLUTIONS 

In order to determine the average properties of unsteady boundary 
layers to order 2 2 , it is evident from equations (27) to (32) that the 
second order (in £q) solutions (denoted by double-zero subscript) must 

be known. These solutions are determined in subsequent sections by ex- 
tending the results of references 1 and 2. 


Differential Equations and Boundary Conditions 

The partial differential equations describing the unsteady flow and 
heat transfer of compressible viscous fluids are given in reference 1 as 

t-pp + - ^x^YY = U '( T ) + Cv »*yxy ( 33 ) 

and 


e T + ~ 


Cv 

PT 9 w 


Cv 

+ (*yy) 


2 


(34a) 
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In reference 2 it was found more convenient to write equation (34a) as 


C v„ 


Cv 0 


+ " ^x®y - pF“ ®rY + £~Te - e w ) ^yy^ ' 

P ^ 


(34* ) 


where 


9-9 



The appropriate "boundary conditions on if are 
* y (X,»,T) = t y (0,Y,T) =* U(T) 

* y (X,0,T) = *(X,0,T) = 0 


(35) 

(36a) 

(36b) 


For the case of an insulated surface as treated in reference 1, the bound- 
ary conditions on 0 are 


0(X,“,T) = 0(0, Y,T) =» 0. 


(37a) 


and 

0 y (X,O,T) = 0 (37b) 

For an isothermal plate, equations (37) are replaced (see ref. 2) by 

8(X,-,T) = ®(0,Y,T) = 0 (38a) 

end also 

0(X,O,T) = 1 (38b) 


Solutions 

The methods of solutions for the insulated and isothermal plate are 
identical to those in references 1 and 2, respectively, but, in each case, 
the method is extended herein to yield second-order results. 

Insulated plate. - For this case the boundary-value problem is de- 
fined by equations “[ 33) , (34a), (36), and (37). For nearly quasi-steady 
flows the stream function given by equation (5) and the temperature func- 
tion given by equation (3a) are substituted into equations (33), (34a), 
(36), and (37) to yield 
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fff(-i5 n ) - 2j ^(0,^) = f(o,5 n ) - 0 (41) 

r ff (0,S n ) * r(~,£ Q ) - 0 (42) 


Introduction of equation (2) makes equations (41) and (42) self- 
consistent (i.e., functions of a and £ n only). For nearly quasi- 

steady flows, £ n « 1 (see refs. 1 or 2); therefore, the functions f and 

r can he expanded as follows: 

f(ff*£ n )3 F(ff) + Zcfoia) + + ■ • • + S^OoM + • • • 

+ £ 0 ^i f oi (CT) + • • • (43) 

r(ff,£ n ) s H(ff) + £ 0 r 0 (cr) + ^r^cr) + . . . + £o r oO^ + • • • 

+ £ 0 ?l r 01 ( ff ) + • • ■ (44) 

Substitution of equations (42) and (44) into equations (39) to (42) and 
collection of terms independent of £ and those multiplied by £q and 

yield the three sets of equations solved in reference 1. Since for 

the present purposes the next higher order terms of equations (43) and 
(44) are needed (see eqs. (27) to (32), e.g.), the terms multiplied by 
o 

are collected and yield 

f oo + Ff bo - ^oo + ^oo = " 2f 0 (2 - f o) + f o( 2ff - 3f 0 ) (45) 
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and 


r 5 0 +Pr(Fr 6 0 


in addition to 


and 


4 *" r oo)- Pr [ 2ffr 6 +2f A I -o - 5S ’ f oo F " f oo -f < f o> 2 ] 

(46) 

f oo(“> - f oo<°) - f oo(°> - 0 ( 47 > 

r oo<°> ■ r oo ( *> = 0 («» 


The function F and its derivatives are available in reference 4, and 
f 0^ r Q and their derivatives are presented in reference 1. Hence, the 

second-order functions f^ and Tqq for the insulated plate are de- 
fined by equations (45) to (48) . 


Isothermal plate . - Equation (45) holds for the isothermal plate as 
well as for the insulated one. However, for the isothermal plate, equa- 
tions (34b) and (38) replace equations (34a) and (37). The temperature 
function given by equation (4) is obtained by letting 


8 - + 2^ U ^: - 6. T 


where 


- B(o) + i^la) + qh^cr) + . . . + C^ 00 (o) + . 


+ ViV a > 


+ • • • 


and 


s(<J,£ n ) = s(o) + t, 0 B 0 (a) + ^ 1 s 1 (cr) + . . . + CqSqqCcO + • ■ ■ 


+ C 5 s (a) 
T1 01 


Substituting equation (4) into equations (34b) and (38) yields those equa- 
tions treated in reference 2 independent of £ and to order and to . 

o {J s. 

To order £q the following equations are obtained: 


fiZRC 



3929 


NACA TN 3886 


13 


* 00 + ^(Fhfo ' = M 2 * h 6 - 3^ - SH*^ + 21^ - SHq) 

(49) 


"So + Pr < Fs oo - 4F ' 0 oo> - Pr [ 2as o - 5S ' f 00 - 3f o s o + 2s o f c. 

- F-fgo 

- I C f o> 2 ] 

(50) 

W> - *W°> " 0 

(51) 

s 00 (co) = 8 00 (0) * 0 

(52) 


The functions on the right-hand sides of equations (49) and (50) are given 
in references 1, 2, and 4. 

Equations (45) to (52) thus define four boundary- value problems for 
the functions £ qq > t qq > ^ qo ’ and s 00’ Tiiese Problems are solved for a 

Prandtl number of 0.72 by a numerical integration method described in 
references 1 and 2. The functions are presented in table I. 


RESULTS FOR AVERAGE PROPERTIES OF UNSTEADY BOUNDARY LAYER 

The following formulas are obtained by substituting the results of 
references 1 and 2, and those of the previous section, into equations (27), 
(28), (29), (30), and (32), respectively: 

C f = (0.6640) ~ (1-306) ffi 2 g~ 5/2 g’ 2 + . . .] (53) 

For the insulated plate case, 


8*= (1.721) / ^^|g“ 1 / 2 + (0.645) m*g 3 / 2 

+ S 2 [(1.986) g - 9/ 2 g ' 2 + (3.420) m*g _5 / 2 g' 2 ] + . . . ^ (54) 

= 1 + (0.8480) ^-[g 2 + (5.051) S 2 g“ 2 g' 2 + . . J (55) 


and for the isothermal plate case. 
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Jc7~x f _____ — t— 

6* * (1.721) <(l + 1.13 S) g" 1 / 2 + (0.1673) m* g 3 / 2 

+ 2 2 [(1.986 + 1.095 #) g-9/ 2 g <2 + (0.3112) m* g" 5 / 2 g' 2 ] + . . A 

(56) 

i . (0.5912) $ _ (0.424) m* 

- 2 2 [(0.8350) $ g _7 / 2 g' 2 + (0.7574) m* g“ 3 / 2 g' 2 +]. . A (57) 

Discussion of Quasi -Steady Term 

The leading terms (i.e., terms independent of 2) of equations (53) 
to (57) reflect the nonlinear dependence of the physical quantities on 
U, in a quasi-steady situation. Under the restriction that g is posi- 
tive, all these leading terms are, of course, positive. 


Obviously, in the completely steady insulated case, when g = 1, all 
the leading terms in the traces equal 1. In the Bteady isothermal case, 
of course, only the factors involving g became unity. In the general 
quasi-steady case, when only the restriction that g = 1 is applied (eq. 

( 21) ) , the magnitudes of the leading terms may he compared with unity (the 
steady case) using a special case of the Schwartz inequality (ref. 5), 
namely. 


For m,n = 0,l/2, the result is 


S 


g2n > g m+n~ 

(58) 

; made, equation (58) yields 


g 2 > 1 

(59a) 

7pz< i 

(59b) 


For m,n = 3/4, l/4, g 3 / 2 g 1 / 2 > g = 1, and therefore, in view of equation 
(59b), yields 


g 5 / 2 > 1 


(59c) 
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If m = n = l/4 , g 1 / 2 g" 1 / 2 > 1, which, together with equation (59b) , 
yields 

g -1 / 2 > 1 (59d) 

If m,n = 5/4, l/4, then g^T 2 g^-7 2 S g^/ 2 ^, therefore, in view of equa- 
tions (59h) and (59c), yields 

g 5 / 2 > 1 (59e) 

Insulated surface . - Inspection of equations (53) to (55), together 
with equations (59), shows that for the insulated surface case, the quasi- 
steady values of C f , 6*, and 0^ are always greater than or equal to the 
corresponding quantity for uniform flow at the average velocity U m . 


The foregoing effects are clarified hy considering the average veloc- 
ity profile. For example, suppose that for half the time, g =* l/2, and 
for the rest of the time, g = 3/2, so that g = 1 (fig. 1(a)); for in- 
compressible flow, the average quasi-steady velocity profile is 



(60) 


The result appears in figure 1(h) ; the dash-dot lines represent the in- 
stantaneous quasi-steady profiles at the low and high velocities, the dash 
line represents the Blasius profile for the average (as if g = l), and 
the solid line is the average profile from equation (60). The velocity 
from equation (60) is greater than the Blasius value near the wall and 
less than that value far from the wall. These profile differences are 
more pronounced if the variation of stream velocity Is more extreme; if 
g « l/2 for nine-tenths of the time and if g = ll/2 for the rest (fig. 
1(c)), then 



The result appears in figure 1(d) . The deviation of the average profile 
from the Blasius shape clearly implies increased skin friction, and also 
suggests an increase of 5*. 

Note that the effects just discussed can he determined without fully 
specifying the function g. All that is needed is a specification of the 
proportion of time during which the various velocity values apply; in ef- 
fect, a probability density distribution for velocity is sufficient. In 
deriving equations ( 60) and ( 61) , the illustrative velocity functions were 
chosen for simplicity; of course, the present analysis would not apply near 
the sudden velocity changes which were postulated. 
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Heat transfer . - If there is heat transfer at constant vail temper- 
ature, the quasi-steady value of 8* from equation (56) is greater than 
the Blasius value unless § is nearly at its maximum negative value, -1, 
which corresponds to a wall temperature of absolute zero. The heat- 

transfer rate q must he regarded as a function of $ and m* (eq. (56)). 
The finer cross-hatched plane of figure 2 denotes q as a function of #■ 

and m* for the average U m (g = l) . Positive q denotes heat flow out 
of the surface into the gas . If g is not always 1, then the first term 

of equation (57), independent of m*, tends to diminish the magnitude of 
quasi-steady heat-transfer rate. The second term, proportional to m* 
and independent of #, always provides less heat flow out of the surface. 
The result is sketched as the coarser cross-hatched plane of figure 2. 

Figure 2 indicates that variable stream velocity results in increased 
heat transfer to the wall only in the case of a cooled (q < o) wall, and 
even then only if m* has a value so that 


m ^ 


(-$) 
0.424 


g 


W* _ n 


(62) 


in cases of negative $ . If # > 0, the condition of zero average heat 
transfer Is 


m 


0,424 


(63) 


Discussion of Frequency-Dependent Terms 

Equation (53) indicates that average skin friction is diminished in 
proportion to the square of the frequency. The dominant factor in this 
effect is the response of the boundary layer to rate of change of accelera- 
tion (fj[(0) in eq. (27)). In respect to skin friction, the quasi-steady 

term and the frequency term are in opposition. Subject to the requirement 
that ffl«l, a flat plate will experience a diminished average drag If 


S2 


2 > g 5 / 2 - I 

(l.306)g" 5 / 2 g' 2 


( 64 ) 


Owing to the appearance of g' 2 in inequality (64), abrupt velocity 
changes would favor drag reduction. 
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Insulated' surface . - The average value of S*, from equation (54), 

2 

will increase with S2 . The adiabatic wall temperature increases with 
S? , which reinforces the effect of the quasi-steady term. 


Heat transfer . - From equation (57), the terms proportional to ffi 2 
all reinforce the effects of the quasi-steady terms with respect to the 

steady value. That is, the <2 ter m in $ reduces the magnitude of q 

as does the quasi-steady term since gV^ < 1, and the negative value of 
o „ 

the 9, term in m* - is consistent with the corresponding quasi-steady 


term since 


?7S 


2 : 1 . 


2 jt 

In the heat-transfer case, the 2 terms of 6 also reinforce the 
corresponding quasi-steady terms, when the groups depending on $ and 

m* are considered separately, Just as in the foregoing discussion of q. 


Results for Harmonic g(f) 

It may often he of interest to evaluate equations (53) to (57) when 
the function g is specialized to have harmonic form 

g = 1 + 6 sin t 

where the restriction 0 < e < 1 is imposed so that flight direction is 
not reversed. 

Arbitrary e . - The various integrals are evaluated for arbitrary 
6 in appendix B in the order of their appearance in equations (53) to 
(57), with the following results: 
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and K and E are the complete elliptic Integrals of the first and sec- 
ond kind, respectively. The results are presented in table II. 
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Small e . - In the event that 6 is very small, to a first approxi- 
mation, each of the foregoing equations for arbitrary 6 may be replaced 
by making the appropriate substitution for IT in the formulas 



( 66 ) 


e Hear unity . - If s - is near unity, then k (eq. (B3)) is also 

near unity, or k' = a/ 1 - k 2 ® 1 - e is nearly zero; and the asymp- 

totic relations 


K(k) * In — ; E(k) * 1 
aJi - s 

apply (ref. 6). Equations (65) become 



Thus, as 6 -* 1, equations (53) to (57) show that the effect of the 

2 

2 terms is greatly accentuated because their coefficients approach in- 
finity from equations (67). For example, equation (53) becomes 

Cj (0.6640) £(1.20) - (0.3918) + * * •] ( 68 ) 

It should be noted that under these circumstances (e -*■ l) the higher order 
terms in S 2 become significant. 
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REMARKS ON HEAT -TRANSFER PROBLEM 

In a previous section it was shown that the heat -transfer rate with 
an oscillating surface differs from that with a surface in steady motion 
(see fig. 2). The relative merits of the oscillations, however, depend 
upon the particular configuration or application. For example, it has 
Been determined that a greater heat-transfer rate to a cooled oscillating 
wall can he obtained for < 0 provided that inequality (62) applies. 

This result suggests that if one wishes to increase the rate of heat ah- 
straction from a gas flow over a cooled surface (as, e.g.. In a heat ex- 
changer), it may he advantageous to oscillate the surface mechanically in to 

its own plane so that the relative stream velocity oscillates. The power 
required to oscillate the plate should properly he assessed to the system, 
because unless this power is less than the additional heat transfer ob- 
tained, that power could he used to increase the energy obtained with the 
surface fixed or in steady motion. Of course, if the oscillations were 
inherent to the system, the assessment would not he required. 

If the power required to oscillate the plate is taken into account 
the comparison is as follows: The excess of power required beyond that 

corresponding to steady flow is given by 3 

J } - c f (u m )] ax (69) 

the wall owing to the oscillation 


q(u m j] ax (70) 


Equation (69) can he evaluated from equation (53), because gC^ is 
given by the same formula as with the exponent of g increased by 

1 wherever It appears in equation (53): 

- 

The formula used for AP accounts for all the power required to 
maintain the plate In flight. If only that power required to oscillate the 
plate in a steady wind-tunnel flow or in a stationary heat exchanger were 

desired, the formula AP » (U - U m ) OX would be used, which is 

less than that given in equation (69) In the amount jf 

Equation (69) is selected for discussion on the grounds that the difference 
cited would usually be charged as an energy loss of a practical airborne 
system. 



and the Increment of heat transfer into 
is 

A<i " -Jo • 

where X is the length of the plate. 
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Equation (70) can "be evaluated directly from equation (57), setting 
g = 1 to define q(U m ) : 


AQ = (0.5912) 


Cp6„pJJ m plx 
Pr 'V Um 


x [(l - g 1 / 2 ) $ + (0.424) m^g 5 / 2 _ +^(S 2 )] (72) 

The question is whether AQ > AP . Therefore, equation (71) is sub- 
tracted from equation (72) to yield 


AQ - AP 


(0.5912) 




* 


+ (0 .424) (l - 2.65 Pr)(g 5 / 2 - l) m* + <3<a 2 )] (73) 

Thus, if AQ - AP > 0, 

~ - 1 > (0.424) (2.65 Pr - l) m* — 1 (74) 

~T7z 

1 - g 1 

If the gas in question is air with a Prandtl number of 0.72, then the con- 
stant is 


(0.424) (2.65 Pr - l) - 0.385 


(74a) 


Thus, for this case, if inequality (74) is to hold, — - 1 must he posi- 
tive. Wow, the greatest value which 0^/S^ can approach consistent with 

heat flow into the surface is the insulated surface value from equation 
(55). Thus, inequality (74) becomes 

0.424 m*g 2 > ^ - 1 > 0.385 m* (75) 

1-572 


or, 

i.ioi g 2 > Big 
1-572 


( 76 ) 
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It is doubtful that any reasonable function of g can be devised to 
satisfy inequality (76) . Inspection of table II shows that inequality 
(76) c ann ot be satisfied by a harmonic function of any amplitude. There- 
fore, for air at normal conditions, the power required to oscillate the 
plate would exceed the extra heat transfer obtained by oscillation. 

However, the constant of equation (74a) is especially sensitive to 
Pr and if, for example. 


Pr< 2^5 = 0 - 3772 ' 

then the constant of equation (74a) would be negative. Thus, in this 
case, when the constant in inequality (74) Is negative, inequality (74) 
can be satisfied for a range of reasonable cooled-wall temperatur es , 
beginning with a minimum value less than 0 W by an amount depending on 

m*. In this circumstance It is advantageous to oscillate the plate, par- 
ticularly for large m*. 

The results of the foregoing discussion may be compared with the 
result of reference 3, wherein it was concluded that oscillation of a 
doubly infinite plate in a fluid otherwise at rest would also result in 
an increased heat transfer to the plate greater than the extra power re- 

e w 

qulred only if — > 1. 

For a stationary configuration, inequality (74) Is computed using 
the AP expression in footnote 3 and is. given by 


e w 

9 


In this case beneficial effects of the oscillation are obtained for some- 
what larger Prandtl numbers (Pr < 0.5 approximately for harmonic oscilla- 
tions) . In any case for Pr -*■ 0 it makes no" difference how AP Is com- 
puted (see eqs. (7l) and (73)), because AQ .will always exceed AP. 
Therefore, for fluids with low Prandtl numbers, such as liquid metals, 
the oscillations will always yield greater heat-transfer rates . 


1 > 0.424 m*[2.65 Pr(g 5 / 2 . - g 5 / 2 ) + 1 - g 5 / 2 ] 

1 - 77 * 


CONCLUSIONS ! * 

The unsteady laminar boundary layer bn a T flat plate in compressible 
flow has been analyzed for the case of time -variable velocity of flight 
with a view to describing the time-average characteristics of such a 
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"boundary layer. Flight velocity is assumed to vary slowly enough so that 
the resulting "boundary- layer flow is nearly, but not quite, quasi-steady. 
The wall temperature has been assumed constant both along the plate and 
in time; further analysis would be required for the case of fluctuating 
surface temperature. 

In order to obtain time averages, the expansions of flow quantities 

= xu* 

U ~ 2 

£ 0 and were taken from previous work, while the terms were ob- 

tained by numerical integration, and are presented in tabular form herein. 

The time-averages of skin friction and "displacement thickness" are 
presented, as well as heat -transfer rate at the surface, and for the spe- 
cial case of an adiabatic wall, the surface temperature. 

A significant amount of information may be obtained without specify- 
ing U(t) beyond the requirements that U(t) remain positive and have an 
average value U m . Each time average involves two groups of terms, to the 

order contemplated in the present report: A time average of quasi-steady 

terms, and terms proportional to the inverse square of the characteristic 
time of the velocity fluctuation (i.e., the square of reduced frequency). 

The quasi-steady terms differ from the values for steady flow at the 
corresponding average velocity, owing to the nonlinear dependence of the 
physical quantities on U(t) . In fact, especially for extreme variations 
of U(t) about the mean value, the average velocity profile is steeper 
near the wall and more gradual in its outer portion than the Blasius pro- 
file which applies at each instant. Thus, in the quasi-steady approxi- 
mation, skin friction, "displacement thickness," and adiabatic wall tem- 
peratures are greater on the average than for the case of constant 
velocity. The magnitude of the part of the heat -transfer rate that is 
independent of Mach "number is less in the quasi-steady approximation, 
whereas the Mach-number dependent part differs in the direction of less 
heat out of the surface into the gas. 

The differences cited are reinforced by the frequency dependent 
averages for adiabatic wall temperature and heat-transfer rate. The ef- 
fects oppose one another in the case of skin friction. 

The various time averages are derived for the special case of harmon- 
ic velocity variation. Large amplitude affects chiefly the frequency- 
dependent terms, greatly accentuating their importance.-' 

The question discussed is whether it would be advantageous to oscil- 
late the surface of a heat exchanger in order to take advantage of the 


must include terms in £q 


(tl, 
2 * 


, - X“U w (t) 

5 1 “ * 


and £q. The terms for 
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increased rate of heat transfer to the wall, and it is concluded that the 
heat-transfer advantage would generally he vitiated by the power require- 
ment for oscillating the surface against the action of skin friction, un- 
less the quantity 2.65 Pr is less than unity. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 5, 1956 
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APPENDIX A 
SYMBOLS 

constant defined by eq. (id.) 
local skin-friction coefficient 

specific beat at constant pressure 

related to stream function for flat plate in steady flow- 

functions related to stream function for unsteady flat-plate flow, 
i = 0, 1, 2 , 00 , .. . 

function related to plate velocity, defined by eq. (20) 

temperature function related to steady isothermal flat-plate flow 

functions related to temperature for unsteady isothermal flat-plate 
flow, i » 0 , 1, 2 , 00 , . . . 

thermal conductivity coefficient 

Mach number 

function related to Mach number, defined by eq. (15) 
constant defined by eq. (3l) 
general exponent 
power 

Prandtl number 

total heat-transfer along plate 
local heat-transfer rate 

function related to temperature for steady insulated flat-plate 
flow 

functions related to temperature for unsteady insulated flat-plate 
flow, i = 0, 1, 2, 00, . . . 

function related to temperature for steady isothermal flat-plate 
flow 
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s,S£ functions related to temperature for unsteady isothermal flat-plate 
flow, i = 0, 1, 2, 00, . . . 

T,t time 

U stream or plate velocity in X-direction (see sketch (t>)) 

U mean velocity 

m 

u relative velocity in X-direction 

u a absolute velocity in x-direction 

X coordinate along surface measured from leading edge 

x coordinate along surface in system fixed in fluid (see sketch (a)) 

Y coordinate defined by eq. (8) 

y coordinate normal to surface 

X ratio of specific heats 

v, 

6 displacement thickness 

6 amplitude of velocity fluctuations 

£ n dimensionless parameter, n = 0, 1, 2, . . . (eq. (2)) 

0 dimensionless temperature difference 

0 temperature 

p absolute viscosity coefficient 

v kinematic viscosity coefficient 

p density 

a dimensionless coordinate, defined by eq. (6) 

y 

a m dimensionless coordinate, tt 

t function related to time by eq. (22) 

local wall shear stress 
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<$ constant defined "by eq. (17 ) 

^ stream function 

2 frequency parameter, defined "by eq. (24) 

to frequency of velocity fluctuations 

Subscripts : 

w evaluation at wall (Y = 0) 

• evaluation in stream (Y -*• •») 

Subscript notation for partial differentiation is used when conven- 
ient. Primes denote ordinary differentiation. 

Superscripts : 

time average as defined in eq. (21) 
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APPENDIX B 


DETERMINATION OF INTEGRALS FOR HARMONIC OSCILLATIONS 

Given that g(x) = 1 + e sin t, the various averages appearing in the 

equations are to "be found. The averages g - ^/ 2 and can quickly he 

expressed as complete elliptic integals hy replacing t hy 2z and tak- 
ing account of the symmetry of the integrands: 

2rt „ At/ 2 


g 


-l/2 * _L 
~ 2n . 


ds 


dz 


vr + e sin t rt 


^ 1 - k 2 sin 2 z 


2K(k) 
it//l + e 

(Bl) 


^T/2=±J^ Ajl + e sin x dx = J ^ ^/TT 


k s sin^z dz 


_ 2/\Jl + e 


E(k) 


(B2) 


where 


k = pe/(l + e) 


(B3) 

and K(k) and E(k) are complete elliptic integrals of the first and sec- 
ond kind, respectively, (eqs. (773.1) and (773.3), and tables (1040) and 
(1041) of ref. 6). 

With help from equations (858.3), (436.00), and (436.03) of reference 
6, these results are obtained: 


g 2 3 -i- 

K 2rt 


Art 

/ (1 + 

Jo 


e sin t) dr = 1 + ~ 


g-V 2 


£ f 2 ” 

“ 2rt 

Jo 


(l + e sin t) - 2 cos 2 ^ dx a -1 + — ===== 

<Vl - e 2 


(B4) 

(B5) 


The remaining averages are of the forms 

Art 


N- 


g 


2 = 1 



NH 


(l + s sin t) 2 dr 


(B6) 
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and also 


N+ i 
2 ,2 
.g S' 


2it 


2 

e 

2* 




(l + e sin x) cos^x dx 


fo 


2ic 


■t 1 - * 2 > h 


N+i 

(l + e sin x) dx 


2it 


N+p- 

(l + e sin x) dx 


2jt 


1 
2 n 


N+2 

(l + e sin x) dx 


(B7) 


Thus, equation (B7) can he' evaluated knowing equation (B6), which may he 
obtained from a recursion formula, as follows: From integrating hy parts. 


f2* 


N-t 


€ i (1 + e sin t) 2 coA dT - , (a ,% 3) 

Jo OO 


rt(2K + 3) 


P2n 


(l + e sin x) sin x dx 


2n 


N+§ 

(l + s sin x) dx 


2n 


N4 


(l + e sin t) dr 


Jo 

Combining equations (B7) and (B8) yields 


(B8) 


I ■xX _ N-S I7+S- / - v N+^- 

(n + |)(1 - s 2 ) g 2 - 2(N + 2) g 2 + (* + I) g 2 = 0 (B9) 
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Substituting equations (Bl) and (B2) into (B9) yields all necessary rela- 
tions for evaluating equations (B6) and (B7): 


7^ = JL JTT7 [4E - (1 - e) k] 


g 5 / 2 - ~ /fl + s [(23 + 9e 2 ) E - 8(1 - s) k] 

- -I Vrrr (, - ^ 

4 


g -7 /^' 2 


15rtVl + 6 (1 - e 2 ) 


(BIO) 

(Bll) 

(B12) 

(B13) 


( - r - H -- E - K ) 


g - 9 /2 gl 2 


2,2 


105*/)! + e (1 - e ) 


[(3 + 16s - 6e 2 ) E - 2(3 + 4s - 3e 2 ) k] 


(B15) 

The foregoing averages are presented in table II as functions of e. 
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TABLE I. - SECOHD-QBDEK SOLUTIONS 


cr 

■P 11 

r 00 

f 6o 

f 00 

r 6o 

r 00 


11 00 

s Ao 

B 00 

0 

1.08350 

0 

0 

0 

1.51298 

- 0.13182 

0 

1.03328 

0 

.1 

1.05201 

.10728 

.00539 

-.45493 

1.48921 

-.13245 

-.01321 

.60292 

.08043 

.2 

.96756 

.20864 

.02126 

-.79778 

1.42573 

-.13329 

-.02649 

.31413 

.12530 

.3 

.84501 

.29952 

.04677 

- 1.04778 

1.33275 

-.13283 

-.03981 

.12345 

.14650 

.4 

.69821 

.37683 

.08071 

- 1.22022 

1.21876 

-.12966 

-.05296 

-.00100 

.15217 

.5 

.53920 

.43876 

.12162 

- 1.32738 

1.09088 

-.12273 

-.06562 

-.08207 

.14773 

.6 

.37789 

.48460 

.16792 

- 1.37922 

.95513 

-.11146 

-.07737 

-.13533 

.13668 

.7 

.22196 

.51452 

.21801 

- 1.38395 

.81660 

-.09577 

-.06776 

-.17062 

.12125 

.8 

.07702 

.52936 

.27032 

- 1.34855 

.67967 

-.07608 

-.09638 

-.19353 

.10297 

.9 

-.05306 

.53042 

.32342 

- 1.27929 

.54802 

-.05328 

-.10286 

-.20673 

.08287 

1.0 

-.16583 

.51932 

.37600 

- 1.18215 

.42474 

-.02857 

-.10696 

-.21111 

.06190 

1.1 

-.25986 

.49788 

.42694 

- 1.06312 

.31232 

-.00329 

-.10855 

-.20687 

.04093 

1.2 

-.33452 

.46800 

.47530 

-.92844 

.21264 

.02116 

-.10764 

-.19421 

.02081 

1.3 

-.38974 

.43162 

.52032 

-.78455 

.12694 

.04352 

-.10439 

-.17382 

.00237 

1.4 

-.42593 

.39069 

.56147 

-.63797 

.05582 

.06272 

-.09905 

-.14709 

-.01373 

1.5 

-.44399 

.34704 

.59837 

-.49498 

-.00077 

.07802 

-.09199 

-.11611 

-.02692 

1.6 

-.44530 

.30245 

.63085 

-.36120 

- .04348 

.08895 

-.08360 

-.08334 

-.03689 

1.7 

-.43189 

.25847 

.65888 

-.24104 

-.07347 

.09541 

-.07435 

-.05113 

-.04361 

1.8 

-.40628 

.21648 

.68261 

-.13766 

- .09225 

.0753 

-.06466 

-.02166 

-.04721 

1.9 

-.37091 

.17753 

.70228 

-.05322 

-.10164 

.09577 

-.06496 

.00333 

-.04809 

2.0 

-.32905 

.14254 

.71825 

.01202 

-.10354 

.09073 

-.04562 

.02295 

-.04674 

2.1 

-.28380 

.11182 

.73092 

.05906 

-.09984 

.08314 

-.03690 

.03695 

-.04369 

2.2 

-.23782 

.08580 

.74078 

.08980 

-.09227 

.07381 

-.02905 

.04565 

-.03953 

2.3 

-.19386 

.06419 

.74823 

.10679 

-.08234 

.06352 

-.02216 

.04973 

-.03472 

2.4 

-.15345 

.04688 

.75374 

.11281 

-.07128 

.05298 

-.'01635 

.05009 

-.02969 

2.5 

-.11804 

.03339 

.75776 

.11065 

-.06005 

.04278 

-.01156 

.04770 

-.02479 

2.6 

-.08819 

.02304 

.76051 

.10287 

-.04934 

.03338 

-.00776 

.04347 

-.02022 

2.7 

-.06387 

.01559 

.76247 

.09172 

- .03959 

. 02507 

-.00486 

-.03828 

-.01613 

2.8 

-.04494 

.01008 

.76369 

.07900 

-.03105 

.01799 

-.00271 

.03263 

-.01258 

2.9 

-.03052 

.00642 

.76452 

.06587 

-.02380 

.01225 

-.00121 

.02698 

-.00961 

3.0 

-.02006 

.00392 

.76507 

.05335 

-.01785 

.00776 

-.00022 

.02171 

-.00717 

3.1 

-.01271 

.00223 

.76528 

.04207 

-.01309 

.00440 

.00038 

.01704 

-.00523 

3.2 

-.00773 

.00134 

.76557 

.03233 

-.00939 

.00200 

.00069 

.01304 

-.00374 

3.3 

-.00448 

.00060 

.76555 

.02424 

-.00657 

.00040 

.00081 

.00976 

-.00260 

3.4 

-.00238 

.00038 

.76566 

.01774 

-.00449 

-.00059 

.00080 

.00714 

-.00177 

3.5 

-.00116 

.00018 

.76573 

.01268 

-.00297 

-.00111 

.00071 

.00509 

-.00116 

3.6 

-.00046 

.00006 

.76562 

.00886 

-.00191 

-.00133 

.00059 

.00355 

-.00073 

3.7 

-.00015 

.00013 

.76579 

.00604 

-.00117 

-.00133 

.00046 

.00243 

-.00044 

3.8 

.00009 

-.00003 

.76561 

.00403 

-.00067 

-.00124 

.00033 

.00164 

-.00022 

3.9 

.00010 

.00014 

.76575 

.00263 

-.00034 

-.00106 

.00022 

.00108 

-.00009 

4.0 

.00019 

.00006 

.76576 

.00169 

-.00013 

-.00087 

.00012 

.00070 

-.00001 

4.1 

.00021 

.00013 

.76565 

.00106 

.00000 

-.00070 

.00003 

.00045 

.00004 



/ foo iff = 2.5490 

r 00 dtr 

= 1.0422 

w 0 

- 0.1633 

/ s 00 
Oq 

= 0.0713 
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TABLE II. - EVALUATION OF INTEGRALS 


e 

*1 

*2 

*3 

*4 

*5 


17 

is 

19 

J 10 

0.2 

1.00748 

.02087 

1.00769 

.02269 

1.020 

.02059 

0.99747 

1.03747 

.02168 

.02050 

.4 

1.03019 

.09673 

1.03299 

.13704 

1.080 

.09107 

.98960 

1.14963 

.11282 

.08680 

.6 

1.06872 

.29227 

1.08539 

.71028 

1.180 

.24999 

.97522 

1.33554 

.43430 

.22032 

.8 

1.12432 

.96153 

1.20016 

7.28120 

1.320 

.66663 

.95129 

1.59356 

2.41368 

.49776 

.9 

1.15938 

2.38511 

1.33236 

63.68454 

1.405 

1.29413- 

.93287 

1.74880 

10.88546 

. 79898 


1.17904 

5.30960 

1.47463 

529.649 

1.45125 

2.2025Z 

. 92025 

1.83273 

45.84304 

1.10872 
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